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SYNOPSIS 


The general problem of ultimate load analysis of 
structures with sections exhibiting bilinear strain hardening 
material property is extremely complex and highly computatire . 
A method of ultimate load analysis of such structures, ideal 
for computer application, is presented. All the conditions 
of plastic analysis are directly satisfied without bringing 
in the concept of collapse mechanism. The problem is 
formulated as a linear programming problem with the load 
factor as the objective function. Equilibrium equations, 
yield criterion and the criterion for rotation compatibility 
are posed as constraint equations. 

The presented algorithm consists of two cycles of 
operation. In the first cycle the structure is analysed 
with the assumption that sections have elastic-plastic 
characteristics. The inelastic rotation and rotation 
capacity of each section are computed and then additional 
constraints for rotation compatibility are imposed. In 
the second cycle, the analysis is done with these additional 
constraints , 

Examples have been included to illustrate the 
advantages of such a procedure. 
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INTRODUCTION 


Most structural materials exhibit considerable 
reserve stren£,th and ductility beyond elastic range. The 
material behaviour in post-elastic region is generally 
nonlinear (Pig. 1,1 a). Rigorous analysis based on nonlinear 
U-0 relaxionship is very complex and time consuming. Even 
the use of a bilinear moment-curvature relationship, which 
IS the simplest mathematical model closely representing 
the inelasxic behaviour of the material; is a highly 
computative process. The development of plastic theory 
based on elastic-perfectly plastic (Pig. 1,1b) idealisation 
of the material behaviour has been only a matter of 
convenience. This assumption of elastic-plastic moment 
curvature relation is justified for certain grades of 
steel. But materials like high carbon steel, aluminium 
and concrete, which have a nonlinear moment curvature 
relationship, can also have very limited ductility. The 
two factors put together may lead to large errors if the 
structures are analysed and designed on the basis of 
assumption of ideal elastic-plastic behaviour v/lth un- 
limited ductility, 

1 .2 Definitions 

According to the elementary plastic theory, as the 
loads on a structure are gradually increased , plastic hinges 
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are gradually fomed at several critical sections. At the 

sections v/here plastic hinges have appeared, the hinges 

continue to rota be under constant moment. The other sections 

continue to take shares of the applied load till sufficient 

number of hinges from causing the structure to collapse by 

1 2 '5 

the formation of a mechanism ’ ’ . This process of transfer 

of moments is called the " moment redistribution " , Complete 

moment redistribution has to take place in order to be able 

to exhaust the post-elastic reserve strength fully. This 

implies that the hinges formed earlier will freely rotate 

to such extenxs as necessary to permit the later hinges to 

form. Such rotations of the hinges will cause large 

curvatures and consequently large strains. This is where 

large ductility plays an important part in plastic analysis , 

Even for frames made of highly ductile materials, it has 
4 ^ 

been recognized ’ that such factors as shear, axial force, 
local buckling, residual stresses and lateral buckling may 
cause the moment curvature relationship to differ markedly 
from the idealised one on which the concept of the collapse 
mechanism is based. In vie\; of the above factors, inelastic 
analysis and design with materials having nonlinear behaviour 
and limited ductility is an important problem, 

1 ,3 Necessity of Present Work 

Several attempts have been made to take into account 
the above factors, piecewise^^"^’ ^ . Nevertheless, the 
existing methods of inelastic analysis and design are subject 
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to the following Imitations: 

1 . There is no conprehensive analytical method to 
evaluate ultimate load providing a direct solution satis- 
fying all the conditions of plastic analysis. 

2. There is no design method adoptable for design 
office taking into account the effect of strain hardening 
and effect of limited ^fuctility. 

3. The methods are largely trial and error 
procedures . 

4-, The existing methods cannot be easily adopted 
to computers, 

5. There is no plastic design method to handle 
cases of partial collapse. 

1 A Object and Scope 

In the light of the above discussions, the object 
of the thesis is the following; 

1 . To propose a method of plastic analysis using 
linear programming where the ultimate load is evaluated 
without the concept of collapse mechanisms » ill the 
conditions of plastic analysis are directly satisfied, 

2. To provide a comprehensive method of analysis 
taking into account the effect of limited ductility and 
strain hardening. An algorithm is presented, 

3. To provide a method for handling the case of 
partial collapse.. In the procedure only a possible 
design solution is obtained as a first step in overcoming 
one of the difficult problem of designer. 
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V/ith the availability of bigger and modern conputcrs 
there is scope for the proposed method to be very effectively 
developed and used as a regular mexhod for inelastic design 
of structures, 

1.5 Proposed Method 

A linear programming approach to obtain a possible 
design solution in cases of partial collapse is presented 
in Chapter V. In Chapter VI, the sane approach makes use 
of the bilinear, strain hardening material property with 
limited ductility to obtain the limit load for plane 
frames, A tv/o step procedure is developed to arrive at 
the ultimate load. 

Material nonlinearity and limited ductility are 
taken into account by formulating constraints on -inelastic 
rotations at critical sections in the structure. In the 
first step of the analysis, an elastic plastic characteristic 
of the section with infinite ductility is assumed. The 
problem is posed as a linear programming problem and the 
ultima be load and the corresponding moment distribution are 
obtained. An expression for the inelastic rotation in the 
above structure is developed and, finally, constraints on 
inelastic rotations are added to the previously formulated 
problem. The ultimate load of the structure is now obtained, 
with these additional constraints. 
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In Appendix I, a flow chart for progranning this 
algorithm is presented. 

1 ,6 Assumptions 

The following are the assumptions made in the 
method proposed. 

1 . The strucxural material may or may not be 
sufficiently ductile to sustain the large local strains 
required for complete redistribution of moments. 

2. Because of the limited ductility and the strain 
hardening behaviour, inelasticity will develop around the 
critical sections where magnitude of moment is more than 
yield moment (M ) , The continuous change of slope in this 
inelastic zone is treated as concenxr'ted at critical 
sections . 

3. The effect of axial force end shear force are 
assumed to be included in the moment -curvature relation. 

4. The structure does not become unstable until 
atleast one section reaches its ultimate curvature (0^) . 

5 . The loads on structure are assumed to increase 
propor ti onately . 

6. Deformations at ultimate load are small so that 
they do not appreciably influence the equilibrium condition 
for the structures. 
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CHAPTER II 


METHODS OP PLASTIC ANALYSIS 

2.1 Introduction 

The elementary theory of plastic analysis or design 
IS based on the satisfaction of three basic conditions. 

1 . Equilibrium 

2. Yield condition 

3. Mechanism condition 

While the first of the above conditions ensures the structure 
to be in equilibrium at all stages of loading upto failure, 
the second ensures that the moment at any section is less than 
or equal to the plastic moment capacity of the section. Under 
the idealisation of elastic perfectly plastic moment curvature 
relation, a section behaves as a perfect hinge under constant 
moment as soon as the section reaches its plastic moment 
capacity. When sufficient number of sections have developed 
plastic hinges rendering the structure an unstable mechaniamj 
the structure is said to have reached the mechanism condition. 

Within the usual assumptions of plastic analysis, an 
analysis or design solution is unique if all the above three 
conditions are rigorously satisfied. The fact that such direct 
solutionsare complex in many cases has led to the following 
theorems^ 

1 . lower Bound Theorem; If an equilibrium distribution of 
stress can be found which balances the applied load and is 



everywhere less than or equal to yield stress, the structure 
will not fail. At most it will have just reached the ulti- 
mate load. 

2, Upper Bound Theorem: For aqf structure, with compatible 
pattern of plastic deformations, the structure cannot be in 
equilibrium if the rate at which the external forces do work 
exceeds the rate of internal \^ork dissipation. 

2.2 Plastic Analysis Methods 

There are several methods existing for plastic 
analysis and design of structures . Of these , the following 
are a few common ones bearing the essential features involved 
in most of them. The essential difference between these 
methods basically lies in the manner of satisfaction of the 
three conditions listed earlier. 

1 , Yield Surface: All possible mechanisms are considered 
and the corresponding limit equilibrium conditions are 
expressed as inequalities in terms of the 'n' independent 
loads involved. Such a set of inequalities define an 
n dimensional surface in an 'n' dimensional load space, 

A state of loading corresponding to any point within this 
yield sui’face defines a safe loading state. The load 
corresponding to a point on yield surface represents a 
state of collapse while the load conditions corresponding 
to points outside the surface cannot be reached, liVhen 'n' 
is two, a geometric interpretation can be given to the 



problem and hence a solution can be obtained geometrically. 
Otherwise one has to resort to mathematical programming, 

Example 2,1: 

Consider the frame shown in Pig, 2.1 (a). There are 
two loads P and H, The plastic moment capacities are 
assumed to be constant throughout the frame and equal both 
in positive and negative directions 

Since the frame in question is a simple one, it is 
easy to see that there are three possible mechanisms. These 
are bhe beam, the panel and the combined mechanisms as shown 
in Pig,2 ,1 (b) , (c ) and (d) respectively. Using the principle 
of virtual work"^ the following inequalities ore obtained 
for the beam mechanism 
for the panel mechanism 
for the combined mechanism 

Each of these limit equilibrium conditions being 
linear is represented by a straight line in Pig.2,1(e). The 
shaded area OABCD defines the safe region for independent 
variation of the two loads. If P/H has a fixed given 
ratio n, as in case of proportional loading, then the 
collapse load is given by point N where P = ;uH intersect 
the boundary of safe region OiBGI). 

2, Method of General Hinge Rotation: Plastic hinges are 
assumed at several critical sections rendering the structure 
a collapse mechanism. The choice of location of such hinges 


P44Mp/L 
P+H 6Mp/L 



1 0 

IS dependent upon the loading condition and also the 
experience of the designer. The following relation is 
established between the number of hinge locations, the 
number of indeterminacies and the number of independent 
rotations; 

m = n-r 

\/here m = number of independent rotations 
n = number of hinge sections 
r = number of indeterminacies rotations. 

The rotations at all the hinges are expressed in teims of 
any of these m independent hinge rotation and the general 
load factor is expressed using virtual work equation4 Those 
independent rotations are given independent variations 
until the least load factor is attained. 

Example 2.2 

The frame shovm in Eig,2,1(a) has three redundancy 
and five critical sections and henco the number of independent 
rotations are 

m = 5-5 = 2 

these two independent rotations are shown in Fig. 2. 2(a). 

The virtual work equation is expressed as 

(H9^ + PGg) L =. Mp(j20^j + j©2“®l| +|2 © 2 I ^ 

p 

for = 0.75, the ultimate load F is given by 

Si _ 2 |Qi| + |e2-^1 I © 2 ! '^I^1+^2| 

Mp ' 0.75 ©1 + ©2 
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Arbitrary values of and 02 a3?e assumed and the corres~ 
ponding value of PL/M is calculated. The result is shown 
graphically in Fig, 2, 2(b), The minimum value of P is 
3.45Mp/L and occurs when Og/Q-i =1 i.e, 02 = ©-j . 

3, Method of Combination of Mechanisms* For complex frames 
with high degree of redundancy, it is not possible to consider 
all the possible mechanisms one by one, primarily because it 
IS difficult t(. cc.'CLive of all the possible mechanisms. 

The method of combining mechanisms is an approximate 

11 12 

procedure presented by Neal and Symonds ’ . In this method, 

a systematic search is made for mechanism which leads to the 
least value of the load factor X , In order to ensure that 
the search at any stage is complete and that there is no 
other mechmiam leading to a smaller value of X , it is 
essential to compute the bending moment at all the critical 
sections to check the yield condition. 

There are, basically three types of mechanisms that 
are considered for plane frames: 

(a) Beam mechanism: In this mode hinges are located only 
in the beams. Pig, 2, 3(a), 

(b) Panel mechanism; In this mode the hinges are located 
only in columns. Pig .2, 3(b) and (o), 

(c) Joint mechanisms; In this mode hinges are located 
only around a joint in all the connecting members. 

Pig. 2. 3(d). 
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In a given structure there are as many basic beam 
mechanisms as there are beams and as many panel mechanisms 
as there are ind<-pendent panels. An independent panel is 
one which can deflect in one unit. Pig. 2, 5(b) and (c) 
show some typical independent panel mechanisms, There are 
as many joint mechanisms as there are joints. A joint is defi- 
ned as onewheietwo or more members me<-t» However^ a certain 
amount of discretion in considering joint mechanisms can be 
used . 


The total number of independent mechanisms 'ra' is 
e^ual to number of equilibrium equations 'e' which is given 
by 

m = e = n-r 

where n is number of critical sections and r is number of 
indeterminacy. 

Example 2,3 

The method of combining mechanism is illustrated by 
means of the frame shown in Pig. 2. 4(a), 

The frame has 6 indetermxnacies . The critical 
sections are numbered as 1 to 14, Therefore number of 
independent mechanisms are equal to eight. QC these 
two are beam mechanisms, Pig, 2,4 (b) and (c); and two are 
panel mechanisms, Pig, 2.4(d) and (e); and four joint 
mechanisms, Pig, 2, 4(f), (g), (h) and (i). The mechanisins 
are combined in a systematic way as shown in Table I, 



T'ABLE Gomb-lnatlon of Mechanisms for Example 
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2,3 imitations of Plastic Analysis Hothods 

In addixion to the methods discussed above there are 

5 

several other plastic analysis and design methods viz. 
method of inequality, method of plastic moment distribution 
etc, The following general observations arc made about all 
these methods 

1 . All the methods arc basically trial and error 
methods. Of the three conditions they have to satisfy 
(viz. equilibrium, yield and mechanism), they uniquely 
satisfy any two and search for a solution among these which 
satisfy the third condition, 

2. The methods of moment distribution and method of 
combination of mechenism are most commonly used. Even in 
those methods after a particular load factor or a particular 
design moment set has been obtained, one has to make a 
complete analysis for the moments at all the critical 
sections; other than the hinge ones, to ensure satisfaction 
of yield condition, 

3. For large frames one has to resort to computers 
in all these methods even though they are not ideally formu- 
lated for computer apolications , 

4. In case of partial collapse the designer is 
faced with the problem of performing another elastic 
analysis of that part of the structure which is still 
statically indeterminate. 
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5 . All these methods are based on an ideal elastic 

perfectly plastic moment curvature relation vath unlimited 

4 5 

ductility. Since this assumption is not alv/ays true ’ , it 
IS more rational to restate the problem with a constraint 
on ductility, 

6. The merit of plastic methods of design basically 

lies in the reserve strength of the structure existing 

between the first yield and collapse. Ignoring the effect 

of strain hardening this reserve strength is normally of the 

order of 15 to 20?^, Since the plastic deformations are 

generally localised, the strain hardening region is reached 

at small deflections. Consequently the section tends to 

develop moments much greater than the idealised plastic 

moment capacity. Pig, 2, 6 shows a typical load deformation 

behaviour of a simple beam test rcjported by Driscoll and 
1 3 

Beedle , In this case a reserve strength of nearly 22 % 
above the theoretical collapse load ignoring strain harden- 
ing IS available. Therefore it is often irrational to 
Ignore the effect of strain hardening which can sometimes 
lead to a reserve strength of 225^ over and above the 15 to 
20?^ reserve strength due to moment redistribution under the 
elastic perfectly plastic idealisation alone. 
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FIGURE 2.2 



TYPES OF MECHANISMS 


figure 2,3 
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(q) Beam Mechonisms 

Frame tor example 2 3 

^ (d) 

Panel Mechanisms 


Frame for Example 2 3 
Figure 2 4 



Joint Mechanisms. 




Collapse Mechanism for Example 2*3 
Figure 2*5 
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CHAPTER III 


PROBLM OP FRAMES V/ITH LIMITED DUCTILITY 

5.1 Introduction 

One of the main reasons for elastic perfectly plastic 
idealisation is that the design procedure bee ernes simple, 

Such idealisations, ho\^ever, lock from one to one correspon- 
dance oetv/een the moments and the curvatures beyond elastic 
limit. There are several materials, whose behaviour can be 
idealised as above and yet limitations have to be specified 
on the ductility. Complications enter when such limitations 
have to be taken into account in design, Nonlinear relation- 
ships arc presently used in Russia and similar trends arc 
shov’/ing in other countries. Design methods based on a 
rigourous satisfaction of even bilinear moment-curvature , 
not to talk about nonlinear, are quite complex. Therefore 
even the bilinear methods are all approximate methods and 
are based on simplifying assumptions. There are several 
procedures proposed to specifically take this fact into 
account. The essential features of these methods lie in 
imposing the ductility limits on the rotation of the plastic 
hinges under gradually increasing loads. This problem is 
commonly referred to as the problem of "Rotational 
Compatibility" , In these design methods the compatibility 
requirement is satisfied by specifying the following 
limitations on the plastic hinge rotations i 
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where 9^ is the inelastic rot-tion at the plastic hinge 
section 1 . This inplics that the inelasticity is concen- 
trated at the cribical sections. 0p^ represents the 
rotation capacity. This iraplies that has to be based 
on a more realistic K-0 relation v/here a one to one 
correspondence between moment and curvature exists, 

3,2 Limitations of the Compatibility Formulation 

There are several limitations in the formulation 
of compatibility condition in the form <4: ®pn.* 
essential to keep these limitations in mind while discus- 
sing the formulation* The limitations are as follows: 

1 * The inelastic rotation 9^ is based on elastic 
perfectly plastic idealisation while 9 is basod on a 
more accurate U-0 relation. Therefore, satisfying the 
compatibility condition as 9^ ^ 0p^ is only an approximate 
and indirect way of incorporating the material behaviour 
into the problem. The extent of this approximation depends 
on the extent of deviation of the actual M-jZ5 model from the 
ideal elastic perfectly plastic one. The closer tho actual 
M-0 is to the ideal elastic perfectly plastic one, the 
leaser the significance of 9 , The farther the actual 

M-0 IS from the ideal elastic-perfeotly plastic one, the 
lesser the significance of 9^. Therefore, the formulation 
IS meaningful only between these extremes, 

2, By definition 6^^ is indirectly the area between 
the idealised M-^ and the actual As the actual 
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gets closer to the ideal one this area reduces and hence 
This IS tending towards more severe compatihility 
restrictions . 

5, The formulation ®pi I'ts meaning 

completely when the actual M-0 is elastic perfectly plastic 
but with limited curvature 0^^ In such cases because of 
the lack of one to one correspondance between the moment 
and the curvature beyond yielding, 9 has no meaning. Only 
one 11-0 model exists for both 0^ and • Yet, compatibility 
cannot be warranted because of limited curvature capacity, 

4. Many of the existing proposals for evaluating 
0p^ are based on parameters derived from the rotations 
observed in tests. These rotations need not always repre- 
sent the real rotation capacity. Consider for example a 
simply supported beam with a concentrated load being tested. 
The actual 11-0 curves of the sections provided in the beam 
IS always bounded between the extremes - an with no 
ductility at all and an yi-0 which is elastic perfectly 
plastic with infinite ductility. In tho test of the above 
beam, at the time of collapse, the maximum rotation that 
can be measured in the case of M-0 with no ductility will 
be zero, In the other case when U-0 is elastic perfectly 
plastic the inelastic rotation that can be measured is 
again zero because the beam is determinate and hence the 
load cannot be increased beyond the level at which the 
moment under the load just reaches the ultimate moment 
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capacity of the section. If a sinalar test is conducted 
on an indeterminate beam the Diaximum rotation that can be 
measured at any hinge section will be equal to the calculated 
inelastic rotation 9^. This again is not necessarily equal 
to the rotation capacity because by assumption the section 
has infinite rotational capacity based on the elastic- 
porfectly plastic with infinite ductility, In reality, 
hov'/evor, the relations are not clastic-perfcctly plastic. 
But based on the above argument, it can be concluded that in 
case of sections whose behaviour is very close to ideal 
elastic -perfectly plastic, large errors may be involved in 
considering the measured rotations as representing the real 
rotation capacities, 

3,3 Inelastic Analysis and Design 

The following are some of the commonly used methods 
for inelastic design of structures of material with limited 
ductility . 

3.3«1 Method of Imposed Rotations 

The method is based on general interpretation of 

the plastic behaviour given by Collometti^^, Since then 

6 9 

it has been established and proposed by Macchi ’ as a 
standard design method for nonlinoar structures. 

The most convenient method for solving inelastic 
problems of the nature discussed earlier is the method of 
imposed rotations. This method assumes inelastic rotations 


as imposed at the critical sections and the statical effect 
of these rotations arc superimposed on elastic effect of 
loads. This method has been practically applied to 
reinforced concrete and prestressed concrete. However, 
the method is quite general and can be applied to any 
material for v/hich bilinear or trilinear E-0 idealisation 
IS valid. The following are the basic steps involved in 
the methods 

1 , To start with, a complete redistribution of 
moments is considered and the load corresponding to this 
moment distribution is calculated. 

2. Por this load, a moment diagram confirming to 
elastic behaviour is determined. 

5, Inelastic rotation at critical sections are 
obtained from known trilinear moment-rotation relationship, 

4, A moment distribution is obtained to cause the 
dislocation in the structure and this distribution is super- 
imposed on the distribution obtained in step 2, 

5. If the final moment diagram thus obtained does not 
violate yield condition, it will be the actual configuration 
at failure. Otherwise another trial with elastic moment 
distribution for a reduced load is required. 

This method of successive approximation is restricted 
to rather simple continuous beams in which the degree of 
indetermiinacy is not very high. 
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3.3.2 Sav/yer's Method 

In this method a bilinear U-0 relation is assumed 
and the rotation capacities are evaluated analytically. A 
bilinear II~0 diagram representing strain hardening is 
assumed to simplify the inelastic analysis and design 

O 

procedures. The bilinear diagram adopted by Sawyer is 
shown in Pig, 7. Por line OA, curvature is given by 


0 - 


M 

W 


(3.1) 


and for any moment between the yield moment and the 
ultimate moment the curvature caused by inelasticity 
IS 




Ml "My 
El 


(3.2) 


where kp is defined as the plasticity factor and is deter- 
mined by the intercept to the left of origin (Pig,3.1). 


Por bilinear analysis the basic problem is the 
evaluation of excess in the curvature and the bending 
angle caused by inelasticity, Por the evaluation of this 
inelastic bending angle, the familiar moment area principle' 
can be used. Based on this fact Sawyer presented a method 
for ultimate load design of concrete structures. The shops 
followed in this method are; 

1 . An elastic analysis is done for all possible 
combinations of loading and envelope moment diagram is 
drawn. 
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2 . These moments arc adjusted to obtain a desirable 
ultimate moment distribution satisfying equilibrium and for 
economy or practicality of designing for such distributions. 

3. The sections are designed and a bilinear moment 
curvature relation based on certain experimental or theore- 
tical rule IS obtained for each section. 

4. for each possible critical combination of ultimate 
loading, using any set of adjusted moments v/hich satisfies 
statics and falls within the M^^'s provided, calculate the 
corresponding inelastic bending angles (k times inelastic 

JT 

moment area) for each plastic region, 

5. Calculate residual moment set by imposing these 
angles to be concentrated at the critical sections, 

6. Superimpose the moments in step 5 on moments of 

step 2, 

7. Check if the yield condition is satisfied at 
all the critical sections. If it is not satisfied, repeat 
steps 4 to 6. 

8. Check the serviceability criterion using elastic 
analysis . 

In brief, the basic justification of bilinear 
analysis and design'’ > lies in taking account 

of all the reserve strength between elastic limit and 
ultimate stage accounting for strain hardening and yet not 
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exceeding the ductility limits. The two methods explained 
'nbove are basically similar, InHacchi's method a trilinear 
M-9 model is used while in Sawyer's, a bilinear ¥i-0 relation 
IS used, Nevertheless, both are based on the principle of 
ensuring ductility limits by imposing only allowable 
rotations as dislocations at the critical sections. 
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CHAPTER IV 
INELASTIC ROTATION 


4.1 Introduction 

As discussed in Chapter III, in problems with limited 

ductility, IX is essential to compute the inelastic rotations 

0^ at the critical sections. This chapter is concerned vnth 

1 9 

the e-valuation of the inelastic rotation , The operations 
are carried out in terms of the residual moments which remain 
at the critical sections when the structure is elastically 
unloaded from the state of collapse to the no-load condition, 

4.2 Residual Moment System 

Consider a structure which is loaded to a loading 
state when n hinges have formed in the structure, and 

It is aboub to collapse. If 

= actual moment at section i at the loading stage Xj^P 

^ei “ pseudo-elastic moment at section i due to load A^P 
(assuming the structure to be behaving perfectly 
elastic during the entire loading history) 

then the residual moment is given by 

Eq,(4,1) referred to hinge sections can be written as 
”l = "ui - *81 

If real hinges are introduced at all the locations of 
plastic hinges, rotation discontinuity occurs at these 
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hinges due to the moment system . These rotations are the 
same as the actual inelastic rotations occurring at the 
plastic hinges under the load 

4.3 Energy of the Residual System 

Consider a member jk between critical sections j and 

k, Pig. 4.1 . If 

m_ = residual moment at end j and 
J 

rajj. = residual moment at end k, 

the .strain energy in the member is given by 

“ak = sil "j “fc 4) (‘*•5) 

where 

= Energy in member 3 k 
L = Unit dimension of length 
^^jk “ ^®Istive stifness of the member 3 k 

then the strain energy for the whole sLructure is 

“r = Sh S + “3 ”k 4) 

where summation is over all the members of the structure. 

4.4 Location of Last Hinge 

In a structure which has developed n plastic hinges, 
Eq.(4.2) gives n known residual moment components. Since 
the residual moments are in equilibrium with aero applied 
load, an equilibrium condition can alv/ays be found relating 
these n residual moments. In general this equilibrium 

/ 
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equation can be written as 

a-jHi^ + b . 2^2 . • ♦ + ^ 1*^1 *** • • • '*' ‘ * * 

+ = 0 ( 4 . 4 ) 

This 'fieans that there are only (n-1 ) independent residual 
moment components while the remaining residual moment compo- 
nents depend on these (n-1) moments. ViTith this the expres- 
sion for rotation is 

.... ( 4 . 5 ) 

In the above expression, represents one of the 
independent residual moments . Therefore for determining ©^ 
values from Eq.(4,5), it is essential to knov/ which of the 
m^ components are independent and which are dependent. Since 
it is required that the inelastic rotation at the last hinge 
section should be zero, the residual moment of this section 
must be a dependent component. Thus, if the location of the 
last hinge is known, all the inelastic rotations at other 
hinge sections can be uniquely determined , 

Consider a structure v/hich has n hinges. lot k 
represent the section where the actual last hinge forms. 

Also assume that all the residual momonts at the hinges act 
in the same direction as the correct inelastic rotations 
(in all the cases, the correct direction for the inelastic 
rotaxions are known as they should be in the direction 
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opposing the ultimate iDoment reaction at the section) so 
tnat the correct 9^ values given by 3Uj,/5in^ are all positive . 
Suppose, on the basic of an arbitrarily assumed last hinge 
at section j , all other rotations are computed and denoted 
by . If the location of the assumed last hinge is wrong, 
the values baaed on the wrong last hinge v/ill be of 
incorrect magnitude with some of them negative in sign. Also 
the rotation 0^ at the correct last hinge k will not be equal 
to Zero. Therefore, ©^ values of rotations must be corrected 
to have the resulting rotations all positive and rotation at 
the last hinge k, zero. This requires a correction -©^ to 
be applied to section k so that 9 ^ - = 0 . Nov/ if -9^ be 

imposed at hinge k, the rotations induced at the other 
hinges can be derived from the principle of virtual work, 

I'rom the equilibrium equation (4.4) the virtual work 
equation is 




+ «n “n = ° 


(4.6) 


As a result, the virtual rotations induced at i due to a 
rotation - 9^ imposed at K is given by 

e- ^ _ g< 

-®k m^ - - aj, 

In general the final correct rotations are 

9 . 


' (Si - si 


(4.7) 
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For this term to be positive or zero for all hinge section 


0 ' 

1 



0 * 




a- 


0 


■k 


(4.8) 


for this condition to be satisfied must be algebrai- 
cally least (l-^rgest negative value), of all values. 


4.5 Summary of Steps to be Followed in Residual Moment Method 

The basic steps involved in determining the inelastic 
rotations by the residual moment method are 

1 , Determine M.,. by multiplying the unit load elastic 
moment oo-efficient by ultimate load ^P, 

2. Find the residual moment nij^ at all the hinge 
locations by subtracting from * 

3. Correct direction is assigned to residual moments, 
so that they act in the expected direction of inelastic 
rotations . 


As a general rule 

adjusted sign = (-).(sign of ) . ( sign of m^ before adjusting) 

4. The equilibrium equations relating the residual 
moraenxs of a mechanism is written. This also determines the 
virtual work co-.0ffici6nt a^^ . 

5. Arbitrarily assume ary section as last hinge. For 
example j and determine derivative 

6. Obtain an expression for 0^ using the Eq.(4.5). 
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7. If the degree of atatical indeterminacy is n, and 
n+1 hinges from (i.e. if the collapse mechanism is complete) 
proceed to step 9 ♦ Otherwise choose as many independent 
cznticol sections as the number of remaining indeterminacy 
(say * • * ) * 


8, Equate inelastic rotation at these sections 

ain^2 

Solve for residual moments '** terns 

of known residual moments. 


9. Substitute all values of in expressions obtained 
in step 7 and obtain values. 

10. If 0^ values are al] positive, then the correct 
rotation is * 8^. If not, find algebraically least value 
of thus determining the actual position of the last 

hinge ' k' . 


11. Then the correct rotations are 
I ®r 

\ st ■ 


Example 4.1 


The method of residual moment is applied to a sinolo 
fixed beam, Eig,4.2(a), to obtain the inelastic rotations at 
collapse . 


The fixed beam is of constant flexural rigidity and 
and constant ultimate moment capacity (+22,0), 
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The load factor at collapse is calculated by 
standard virtual work methbd 

\ _ (1 + 3/2 + ^) X 22 

Pig,4.2(b) represents the pioraent diagram at collapse and the 
elastic moment ^ for load factor = 3.3 is shown in 

cl 'U 

Pig. 4. 2(c). The residual moment calculation is given in 
Table 4.1 below 


TA'RLP 4.1 Residual Foment Calculation for 
Exen^ple 4.1 


Section 

1 

2 

3 

Plastic Moment 

-22.0 

22.0 

-22.0 

Elastic Moment 

c 1 

-29.3 

19.57 

-14.65 

V - Mei 

7.3 

2.43 

- 7,35 

Residual moment with 

adjusted sign m^ 

7.3 

-2.43 

- 7.35 


Evaluation of 6 | values ; 

With reference to Pig . 4 . 2(e), the equilibrium equation 
can be expressed in terras of m^ as 

2m^ + 3^2 + HI 3 = (4.9) 

The strain energy is given by 



Assuming section 1 to be location of last hinge > such that 
m 2 and m 3 are independent while ra-j is dependent^ Eq.(4,9) 


can be T-/ritten as 


rai ~ ®2 "* 


Therefore 


- - 3 . J!l - ^ 1 

— “*■ a ^ "TT 


and 02 = 


isf " "SIT ■* ^ 3 ^ 


and ©2 = 


Ell 


( 4 . 10 ) 


Substituting for and from Tabic 4.1 


* 0,0 


ft' > 10 ♦97 L 

®2 EeI^ 


(4.11) 


®3 = - 


11 .0 I 


Since the assumption that last hinge forms at section 1 , 
yields negati'VB inelastic rotation at section 2 and 3; 
the assumed location of last hinge is incorrect. 


10.97 I 


11 L 

nr 


hence the correct location of last hinge la section 3 where 
Q^/a^ IS minimum. The correct inelastic rotations using 
Eq. (6,7) are 


©^ = 0 + 


= - 


31113 L 22.05 


( 4 . 12 ) 


^ ■’^”^3 ) 


22,20 

EX 


©3 = 0, 


(4.13) 




rr': 




' - 

j stiffness = -t- K 
El 


If 

t^U. €,^jU-C(^ 




JK 


(ci) ^ Typical Member withHesidual Moment 
Figure 4 1 
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(c) - Met 


Elastic Moment Due (d) Collapse Mechanism 

to 20 A y 



(e) Direction of Inelastic Rotation 


Figure 4 2 Fixed Beam for Example, ,4-1 . 
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OHAPTEE / 

APPLICATION OP LIESAfl PEOGRAI'MIHG TO PLASTIC 
ANALYSIS ANL DESIGN 


5.1 Introduction 

The use of linear programning approach to plastic 

24 

analysis was first pointed out by Charnes and Greenberg , 

25 

Later, Dorn and Greenberg developed systematic procedure 
and illustrated it through a simple indeterminate truss 
example. Many other investigators^^’^'^ 
discussed the relationship between limit analysis and 
linear programming. The main advantage in plastic analysis 
of structures by linear programming is that highly indeter- 
minate structures can be analysed easily by a digital 
computer , 

5.2 General Problem of Linear Programming 

A linear programming problem can be mathematically 
stated as 

Minimise or Maximise Z = + •«« + (5.1) 

Subject 10 constraints 
n 

Xa. TX (^or^or=:)b ; i = 1,...m (5*2) 

and 0 0 = l,..,n (5.5) 

where , J = 1,.,, n are variables of the system and 
(5,2) and (5.3) represent constraints on interaction of 
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system parameters X . The non-negativity constraint 

J 

expressed by Eq.(5,3) cah be handled in an easy manner in 
simplex method and hence has not been included in C5»2), 

5*3 Standard Form for Linear Programming Problem 

By addition of certain variables (called slack and 

surplus variables) and some suitable transformation for 

those X which are unrestricted, the above linear programing 

problem is reduced to the form 

n 

minimise Z == X (5.4) 

D=1 ^ ^ 

Subject to the coiistraints 
nv 

X a. _ X = b 1 = 1,2.,.m (5,5) 

j—-] -'-J J -‘• 

where nv is the total number of variables = 1 ^ 0 , of design 
variables + no, of slack variables -i no, of surplus variables 
and X^ ^ 0 3 = 1 , . , , , nv (5.6) 

5,4 Phe Simplex Method 

21 22 

The simplex method ’ is a two phase procedure for 
finding out an optimal solution of linear programming problems , 
A linear programming problem may have 

(a) no feasible solution 

(b) an unbounded solution 

(c) finite number of optimal solutions 

(d) infinite number of optimal solutions. 
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Phase I provides an initial basic feasible solution 
to the sjstern of equations expressed by (5.5) and ( 5 . 6 ). If 
the constraint equations are inconsistent, the problem has 
no feasible solution. This coiidition is reflected at the 
end of Phase I. 


Phase ir uses the basic feasible solution obtained 
in Phase I as a starting point and finds either a optimum 
solution or yields the information that the design space 
is unbounded. 


5.5 Application fo Plastic A.nalysis 

A general plastic analysis problem can be written as 
Maxinise ; Z = X (5.8) 

Subject to the following constraints 
(a) Equilibrium equation, 


®i,n +1 ^ ^ 

1 = 1 , , . . e 

where e = number of equilibrium equations 
and n = number of critical sections. 


(5.9) 


(b) Yield Condition 

6 4 


3 = 1 , . . . .n 


( 5 . 10 ) 


where ' s 

X is 


are the moments at the critical sections 5 
load factor to be maximised; 
s are the constants which are obtained from 
statics f 
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and and M are the po&itive and negative plastic moment 

P J P J 

capacities respectively of critical section . 

Using the transformations 


~ ^ ^n+1 

(5.11) 

the Eq.(5.9) and Eq,(5,10) are transformed to 


n+1 

5r a X = B 
j=1 J ^ 

i ~ 1j«.i e 

(5.12) 

and 


d = 1 , . . . n 

(5.13) 

n 

where = X) a, ^ 

1 j=1 


The plastic analysis problem in standard 

form IS 

Minimise ; Z = - 

(5.14) 

Subject to the constraints 


n+1 

H a, , X. = B. 

' ' i = 

(5.15) 

x, + 

(5.16) 


J — 1 ? • • 

> 0 d = 1 2n+1 (5.17) 

where arc slack variables added to inequalities in 

Eg. (5. 13). 

5.5.1 Nature of Solution 

There are (n+e) constraints expressed by Eq,(5,15) and 
Eq,(5.16) and 2n+1 variables in the system. Thus atmost only 



42 


(n+e) of these 2n+l variables v/ill be non zero in any basic 
feasible solution. Again, since in physical problems, the 
load factor ^n-.1 cannot be zero, it \/ill be a basic variable 
in the optimum solution, For other variables the following 
condition can occur; 


(a) Some are in basis but the corresponding slack 

variablos are not in basis, then = 0 

and from 2q. (5-tl6), 

X, = - M" 

D PI PO 

Using Eq.(b.11 ), moment at section 3 is 


“3 = V 


Hence, bhc moment at section 3 will be equal to bhc positive 
moment capacity of section. 

(b) Some X are not in basis, then the corresponding slack 
J 

variable \;ill be in basis i.e, 

Xj = 0 ana 

again using Eq,( 5 . 1 l), 



Hence, iii this case moment at section j equals negative 
moment capacity of section, 

(c) Both X^ and are basic variables and both arc 
not zero. In such a case moment at critical section will 
be less than the plastic moment capacity of the section 3 , 
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(d) Both ond are in hasis but one is zero, i.e., 

the basic feasible solution is degenerate. In such a case 
tno moincnt at section 3 v/ill be equal to the plastic moment 
capacity , 


Assuming that there are K number of sections for 

which cither X-, or X , . , is not in basis, the remaining 

j n+j+i 

2(n~K)+1 variables will all bo in basis. However, the total 
number of varir'^bles in the basis must be equal to the number 
of constraints and hence 


2 (n—K) +1 + K = n + e 

orK = n- e + 1 ( 5 . 18 ) 

Bvidently (n~e) equals the number of indc-tcrminacics of the 
structure considered. Hence, the following stoteraent is 
made regarding the nature of a optimum basic feasible 
solution of a plastic analysis problem. 

"Any optimum basic feasible solution will alv/ays have 
the moment equal to the moment capacity of critical sections 
atleast at (n-e +1 ) sections. In other words if ra is a number 
of indeterminacy, simplex algorithm v/ill always show formation 
of at least (m+ 1 ) hinges in optimmi solution. However, all 
these hinges may not be involved in the collapse mechanism 
(as in the case of a partial collapse)" , 

As it 18 easily seen, if condition (d) exists, number 
of hinging section will be more than (m+ 1 ). 
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Hus property of optinun solution can be used to 
obt.an a possible design solution in cases of partial 
e ell apse . 

Ibconple 5.1 

Consider the sin^e storey portal frame shown in 
Pig, 5. 1(a), The frame has five criticrl sections. All 
members are of uniform and equal plastic moment capacity , 

The equilibrium equations are expressed as 

-Mj + 211 = 1 50 A 

-M^ + M2 - = 150A 

and the yield condition is 

-100 4,11^^100 (5.18) 
1 = 1 , , . .5 

Substitutiiig - 100 

and X = Xg ( 5 . 19 ) 

the plastic analysis problem is expressed as 

Maximise : Z = Xg 
Subject to the constraints 

X2 - 2 X 3 + + 150 Xg =0 

+ X4-X5+150Xg = 0 

X^ 4200 i =1 , . .5 

X ^>0 1 = 1,.,, 5 


X, -Kg 


(5.20) 



adding slack variables, the problem in standard form is 
expressed as 
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Minimise Z = -Xg 
Subject to 


X2-2X3+X4+150Xg 
X^-X2 +X4_X5+150Xg 


X. 


X. 


X3 


X. 


+x 


7 


+X 


8 

+X, 


+X 


10 


^5 


+X 


11 


= 0 
= 0 
= 200 
= 200 
= 200 
= 200 
= 200 


( 5 . 21 ) 


v/here X^, Xg, Xg, X^q and X^ ^ are slack variables. 

The solution of this linear programming problem obtained 
by simplex algorithm is 


X. = 0 ; X^ = 100: X^ =: 200 

’ ' ^ ^ (5.22) 

X4 = 0 5 X3 = 200 5 Xg = 2 

by back bransformation using Eq.(5,19) 

M.| = -100 ; Mg = 0,0 ; M3 = 100.0 

= -100.0} M5 = 100.0 (5.23) 

Load factor =2,0 


This moment distribution is shov/n in Pig. 5. 1(b) 


46 


5,6 Applicption to Plcstic Design Problems 

A plastic design problem can be stated to be as 
"Given the ultimate load, the structure has to 
withst'^nd, to find out the plastic moment capacity of 
vrriouo nembers in such a way that cost or weight of 
structure is minirium" , 

The forraulation of plastic design problem is illustrated 
by tht; following example. 

Example 5 .2 

The frame shown in Pig, 5. 2(a) consists of three 
members. The weight of members are assumed to be propor- 
tional to lbs plastic moment capacity. Further, to simplify 
the procedure the moment capacity cf two sections are as.^umed 
to be equal. With these assumptions the objective function 
for design problem can bo writtv-n as 

Minimise Z = 4Mp-| + 6 Mp 2 (5,24) 

Constr 'lints on Mp^ and Mp 2 J 

The possible critical raeclionisms for problem are 
shown in Pig, 5 , 2 (b) , (c ) , (d) , (e) , (f ) and (g). The limit 
equilibrium condition for these mechanisms are 


for mechanism 

(b) 

4Mp2 > 3 

for mechanism 

(c) 

4Mp, > 6 

for mechanism (d) 

2Mpi+4Mp2 ^ 9 

for mechanism 

(e) 

2 Mp^+2Mp2 > ^ 

for mechanism 

(f) 

2 Mp^+2Mp2 ^ 3 

for mechanism 

(g) 

4Mpi+2Mp2 ^ 9 


« « « 


(5,25) 
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Out of those constraints if the first four r'-.lations nrc 
satisfied, the rest tv/o \all obviously be satisfied. Hence, 
finally the raininiurn v/eight deeign con be expressed in 
linear programming form as 

Minimise Z = 4M-j + 6 I.I 2 
Subject to the constraints 


4Mp 

>6 


4Mp2 

>3 


2 Mpl + 4Mp2 

>9 


2 Mpi + 2 np 2 

^6 

(5.26) 


The solution of this problem using simplex algorithm is 

Mpi = I ; V " i‘ ' ^ 

5.7 "Problen of Partial Collapse 

Por more complex design problems, it is very 
difficult to conceive of all the possible mechanisms of a 
plane frame. Due to this fact the designer is forced to 
select some preliminary moment capacity of members and 
adjust then to obtain a desired load factor. However j 
if the preliminary plastic moment capacities lead to a 
partial collapse of the frame, then it is a very difficult 
task to obtain a bending moment distribution satisfying 
yield condition by classical methods of plastic analysis. 
This difficulty can be easily overcome by using linear 
programming formulation of analysis problem. Thus linear 
programeung can be very conveniently used tp obtain a 


I 
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better possible design. The sxeps suggested in this 
direction arc. 

1 , Obtain moment capacities of beans based on 
bean ncchanisns and of colurrns based on panel mechanisns . 

2 . Using prelininary rionent capacities obtained in 
step 1 find out the collapse load factor 'A' and the 
corresponding Eioment distribution . If the 'Aa' 
actual lo^d factor rcquirt,d for the frame, then one 
possible design solution for frame is 

M,,_ = ^ X Maximum moment in Member 3 
pj 

where M is plastic moment capacity of member j. 

J 

It can be easily seen that this v/ill lead to a 
better design than obt ^ined in stop 1 . 
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GHAPT'H VI 

PROPOSED LDIIT THEORY ACCOUNTING 5TRAIH 
H />KDENrjG ARP DUCTILITY EPFEOT 

6.1 Introduction 

.-3 mentioned in earlier chapters, limit analysis 
b:ised on bilinear relation with strain hardening effect 
IS quite complex. In this chapter, a method of limit 
analysis using linear prograiiiming technique and accounting 
for strain hardening effects with limited ductility is 
proposed. It is an automatic procedure v/hcre in the three 
o 'iiditions, namely 

(a) equilibrium (b) yield condition and (o) rotation 
c jHpatibility are sinultancously sitisfiod. The method is 
ideally suit^^d for computers. 

6.2 Moment-Curvature Relationship 

In order to define a bilinoar relationship 
(Pig,6 . 1 (a) ) , four quantities, El, CiC and B are quire d. 
For ihe t-'pical bilinear curve in Pig, 6, 1(a) line OA, 
through origin, represents the elastic line for which the 
curvature is given by 

0 = ^ ( 6 , 1 ) 

vfhere M is the moment at the point under consideration, 

Por any moment M-j , between and 11^., the total curvature 



51 


I ! ! K/,? UU, „ 

Turt II,,,, ,^y 



c.'j.i be consndcrGd to consist of t\/o componv-nts 

(a) pseudo-elastic componenx, M^/EI 

(b) excess in curvature, 0^, caused by inelasticity. 

Erom Pig, 6, 1(a), the V'-'luo cf inelastic curvatur- 0^ can 
be expressed as 

Tb - M 

’ib^re B is called the ductilixj^ factor -^nd is given by 

0 

B = ^ (6.3) 

and C5<13 the ratio of moments and My corresponding to 
point Pig, 6,1 (a) . 

6,3 dotation Capacity *0 ' of Critical Sections 

By definition, rotation copicity 6 is a mu-asure of 
bhc ability of a section bo uithstand an applied inel ^stic 
rot-^tion. Hence a limitation on 0^^ also implies a Imita- 
tion on ductility. As the loads on blic structure are 
gradually increased after a critical section reaches its My 
value, inelasticity gradually spreads around the critical 
section. At any instant the length of this inelastic region 
IS referred to as "inelastic length, 1", Pig, 6, 2. 

Pig, 6, 3(a) and (b) shov/ a typical moment distribution 
around a critical section and the corresponding curvature 
distribution. The inelastic region is indicated by the 
shaded region, The rotation capacity 0p of a section can 
be expressed ss 
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0 . 


P 



(6.4) 


v/hure 0 ^ rcprcoents tho inGl'i.'ijtic curvatare distribution 
within the length 1. The quantity 0 corrt, spends to the 
shaded area in Tig, 6. 3(b). UvSing Eq.(6.3), bhe cxprosoion 
for 9^ c 'n be \/ritxon as 


9, 


o 


( ^i:^ ) ( — wr-^) ^ 


El 


assuning the section properties, , B, El and M to be 

i/ 

the sai'ic throughout the inelastic length 1, 





1 

IT 



My) dX 


(6.5) 


The dopiinant p'^raneters influencing rut-^tion 
c ipncity rare 

(a) hinge length, 1. 

(b; the roLio of ultiri'^to curvaturi- 0 ^ and yield 
curvature 0 , B, 

w 

(c) elastic stiffness, El 

(d) shape of the moment diagran around the critical 
sections ♦ 

(e) bhe ratio of moments and Ily, • 

While o(, B, El and depend on the properties of 
the sections, hixigo length 1 and morjcnt distribution depends 
on the loading pattern and behaviour of the structure as a 
whole. Hov/ever, the parameter El cancels from both sides 
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of the coriTpatibility condition 0 ^ ^ 9p^ and h.-rLce its 
absolute value liae no influence on the problctn of conpa- 
tibility . 


6,4 Linit Analysis of Plane Prancs 

Basefj on the assunptions seated in Chapter I and the 
bili„c::r i'~0 curve (Tig, 6,1) for critical sections, the 
projlati of lirnit nialysis is s-cated as 

" G-ivcn the bilinear U-0 property of the critical 
secticas, to fiid the maxinuu lo d which the structure can 
./ith&c"'nd without overstraining the naterial" , 

Mathematically, this problem can be sbitcd as 

Maximise • Z = Load factor A 


Subject to the constraints 

(a) Eopulibriun eouations expressed os Eq,(5.9) 

(b) Yield criterion (Eq.5,IO) 

(c) Eot-tion compatibility specifying that 


0 


1 



(6,9) 



= i 


1 




P 


where ipi 2 ,**^p dvsnote the sections where inelasticity 
has been at'coined. 


for an analysis problem, the configuration of 
structure j loading pattern, and section properties such 
as El, 04 and B are assumed known for all critioal 
sections. Hence, all the constraints discussed earlier 
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can be expressed explicitly in tcrne of raoncnts ct 
critical sections; load factor A and other parameters. 

The conplGxity of the problem lies in the formulation of 
constraints in Eq.(6,9) where explicit expressions for 
0^ and 0^^^ cannot be obtained unless the load factor and 
the corresponding noment distributions arc knov/n. This 
necessitates the adoption of sorae iterative procedure for 
analysis . 

The proposed method consists of finding the ultimate 

load and corresponding moment distribution, first without 

considering the rotation compatibility constraint. Next the 

1 9 

iiethod of residual moment is used to obtain an expression 
for inelastic rotation at critical sections. The proper 
sign convention, as explained in Chr pter IV, assures that 
oil the inelastic rotations will always be positive. Having 
obtained the expression for inelastic rotation, the rotation 
capacities at critical s-^ctions arc calculated. The rotation 
compatibility constraint is now introduced in the linear 
progromning problem already formulated. In the second cycle 
of operation, the ultimate load is found for this new 
nathenatical model of ultimate load analysis problem. If 
more accurccy is desired, a third cycle can be carried out 
v\'ith slight modification in 0 values of the rotation 
compatibility constraint only. 



6.5 Sumiiary of Steps in the Proposed Analysis Procedure: 

The basic stt.ps involved in determining the ultimate 
load for inelastic strucxurc-s are: 

1 . Perm tnc equilibriun equation and the constraints 
on I'l Clients 

2. 3y suitable transfornation (IIq.5.11) and addition 
of slack and surplus v'^riables, express the constraints in 
the standard forn (Eq.5.7). 

3. Maximise, using linccr programming, the load 
factor Xy' subject to constrvaints formulated in step 2. 

4. Obtain, by any elastic analysis procedure, an 
elastic nornent distribution for an unit lo'^d on the struc- 
ture . 

r 

5. Calculate pseudo-elastic moment at load ^^P, 
where P is unit load on structure (assume th- structure to 
beh'^VG in a perfectly elastic manner during entire hisxory 
of loading). 


6. Using the melhod of residual moment distribution, 
obtain expressions for inelastic rotation at those critical 
sections v/hcrc ultimate moment has developed during step 3. 
Thus 


e. 


P 

= r 

j=i 




^i,n+1 ^ 


( 6 . 10 ) 


where j's are only those sections v/herc inelastic rotation 
IS not zero and ore constants. Those expressions 

for inelastic rotations are p in number, 



v/here p = I'lur^ber of inde.terriinacy . 

7. Using the Bq,(6,3) and Eq.(6,4) and raonent distri- 
ooxion obtained in step 3, obtain the rotation capacity of 
the critical section. 

8. Append the consbraint ^ 9p^ to the already 
formulated constraints. 

With the sign convention used in step 6, both 
and Gp^ should always be positive for consistency in de- 
formation, However, if rotation compatibility constraint 
is 9^ 0p^ and load factor is maximised the final moment 

distribution can give a negative value of 0^. This in- 
consistency in deformation is removed by introduction of 
additional constraints expressed as 

G^ > 0 (6.11) 

Thus the irathema bic^'l model for limit an''lysis problem is 
Maximise Z = load factor X 
Subject to 

(a) Equilibrium equabion 

(b) Yield criterion 

(c) Rotation compatibility, Eq,6.10 

(d) Constr'ant ensuring consistt_ncy in rotation, Eq, 

( 6 . 11 ) . 

9. Obtain the maxinum load factor and the corres- 
ponding nomunt distribution using simplex algorithri. 
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10, If more accuracy required, the 0^^^ values in 
constrnnt Sq,(6.9) arc changed corresponding to the moment 
distribution obtained in step 9, and the inequality sign in 
constr'^int (6,9) for these i 'vhere the nonent has become 
Iocs than is changed to equality. The corresponding 

expression in Eq,(6.1l) is rciioved .and simplex algorithm is 
-g:in used for finding the ultiriate load with respeeb to 
ncv'/ly formed coiistr''ints , 

Examples 

The proposed method of ultimate load analysis is 
illustrated by three problems 

(a) a fixod beam 

(b) single storey, single ba,-’' porto.l franc 

(c) single storey, t\/o bay portal fr'^mc , 

The moment curvature relations for the critical 
sections are assumed to be same in positive and ncgLitive 
moment. However, thw method is g^.ncral and can be applied 
to V 'nation in this property, 

Exomiolc 6,1: 

It IS required to find out the ultimate load and 
corresponding moment distribution for a fixed end bean 
shown in Fig, 4, 2(a). The bilinoar M-^ property of critical 
sections ore defined as 

Mu = 22 

B 6,1 


o<= 1 .1 

El = 2900 . 
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IiiL-l &tic Curvature! 


Por any moment I-l, between II^. and 

inelastic 

curv'^turc is 


P ' M - 3L, 

gj ~ (£_riL\ ( y) 

'' oi _ 1 ^ El ’ 


50(M 20) 

El 

( 6 . 12 ) 

First Cycle of Analysis • 


The problem of uitiinatc load analysis 

IS 

Maximise : Z = A 

(6.13) 

Subject to 


(a) Equilibrium equation 


2M-,-3M2+M3+40X= 0 

(6.14) 

(b) Yield Criterion 


^22 4 \ 422 

(6.15) 

1 = 1 , . . . 3 


The optimum solution of this problem is 


- - 22 , Mg = 22 , M 3 = - 22 , A= 5.3 

(6.16) 


This rioraent distrihution is shown in Pig .4.2(b). 
Prom Table 4.1, it is seen that 

rii = + 8, 88 A 

EI2 ” ~^ii^ "f" 5 • 92 A 

+ 4 .44 A 


59 


Subs'Ll tuting these vuluos in Eq,(4.12) and Eq.(4.13) for 

iiielrstic rotations 
"jM, - 13.33 X 
01 „ ^ — 

-Ml - 6 M 3 - 53.28A (6.17) 

©2 _ _ 


Ro b .ti on Cap-Qi ty; 


j X 0.046 3 

0,0 = ^ y- 0.138 :> 

P'^ 

0p3 = ^ X 0.092 : 



2 

g^.jO 

50 ^ 

r\ 

6.90 

TUT ^ 

1 

* 

50 ^ 

0 

4.60 

EX ^ 

2 

= El 


( 6 , 1 s) 


Rotation GomiJfltibility Constrrint’ 

Prom Eq.(6.17) -ind Eq.(6,18), rot'''tion compotibility 
consxrsint -^re 


-3M, - I 3 . 33 A < 4.60 

^ V , (6.19) 

-3Mi - 6M3 - 53,28 ^ 13.8 

and for 0 Lo be positive 


“ 3 II 3 - 15 . 33 A^ 0 
-. 3 M 1 - 6 M 3 - 53.28A^ 0 


(6.20) 


So c :> nd Cycle of Analysis ; 

The maximum value of A, subject to constraints 
oxj^rcssod byEq.(6,14), (6,15), (6.19) and (6,20) is 
tabulated in Table 6.1. 
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T-flBLE 6,1 


Section 

1 

2 

3 


Moment, M 

-22 

18.10 

13. 

,95 

Inelastic 
Rotation 0^ 

4.58 

0.0 

0, 

.0 

Rotation 

Cap'- city ©p^ 

5.04 

0.0 

0 

,0 

Ultimate load 

factor / 

= 2.8 




The rbovc solution hos been obtained using a standard, 

' slisrc ' subrouTinc for nininizaticn of a function by 
sinplcx algorithm . It is seen thot the actual load 
factor IS 15/^ less than that predicted by mechanism 
method of analysis. 

Example 6.2 

It is required to find out the ultimate load for the 
frame in Eig,6,4(a). The Section properties are 

= 22 

B = 3.6 
0 < = 1 .1 
El = 2900 

Inelastic Qurvaturv:; 

For any moment M between M and M, , the inelastic 

y ^ 

curvature is 


= S(M-20) 


( 6 . 21 ) 
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First Cycle of Analysis 

The ultinate load analysis problen v/ithout 
conpatibility constraint is 
Maximize : Z = /\ 

Subject to the constraints 

Mg - 2M5 + + 20 A = 0 

- Mg + ~ + 30 A = 0 

-22 < 4 22 

1 = 1 , ... 5 ( 6 . 22 ) 

The solution of this problem is 

= -22. O5 Mg = 13.2; = 22.0; 

M 4 = - 22 . 0 ; M 5 = 22 . 0 ; A - 2.64 ( 6 . 23 ) 

Inelastic Rotations 

The distribution of and are tabulated 

below in Table 6 . 2 . 

TABIiF 6 . 2 Residual Moment Calculation 


Section 

1 

3 

4 

5 


-22.0 

22.0 

-22.0 

0 

• 

CVJ 

CM 

“ei 

-18.2 

17.63 

- 25.7 

26.95 


- 3.8 

4.37 

3.7 

-4.95 

Residual moment 
n with adjust^ 
ea sign 

-3.8 

-4.37 

3.7 

4.95 
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The positive dircctioii of residual nODents rn^ arc shown in 
Pig. 6, 4(d), The e nailibriuri epuations in terms of residual 
. .’.•'ents are 


h-] + ^2 + = 0 

+2n^ + 2n^ + = 0 (6.24) 

Pron Eq. (6 .24 ) 

= -4.85 (6.25) 

^-ssuning the last hinge to form at Section 3, the inelastic 
rot-tions ab Sections 2 and 5 should be zero and hence Ei2 
and are the dependent coraponents of residual lajiaent 
vS stea, 


Thus nsnie, Eq,(6.24), 


9 nip 0n2 



^ = -I 

01-15 


(6.26) 


Rotation calculation for this problen is given in Table 6,3. 

Expressions for residual nonenbs in torus of ectual 
momem: at sections 1, 3,4 and 5, where hinges have formed 
are 

+ 6.9 'X 

= -(H^ + + 26.859) 

+ 9.75 A 

EI5 := -.(M3 - 10.2 A) 


(6.27) 
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TABLE 6,3 Rotation Calculation of Exanple 6,2 


Mcnber 

U . 

ID 

1 auj, 

®i°as7 

e' -SJIe 

4 

e' 

5 -9115 

1 - 2 

2 2 
ri-j +ni'|n2+ia2 

-n2 =1.15 

-ni.| -2^2=1 3 . 50 

-2m2-1 3.50 

2-4 

2 

n2+Ei2n^+in^ 

-2m2~ni^=6 ,0 

-m2+n^ =8,55 

-202-^4= 6.0 

4 - 5 

2 2 
r^4+'^i4^5+°5 


2n^+ra^ =11 .35 

n4+2Ei^=13.60 

Total e 

1 

7.15 

34.4 

31 .10 
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Using the'sc expressions for residual nonents, inelastic 
rotations are; 

9 ^ = + 2M4 - 3M5 + 77 .70 A 

0^ = 2M^ + 6M^ - 4M5 + 1 1 3 J A 

83 = 3M^ + 4M^ - 6M5 + 120.9 A (6.28) 

The 9 values are all rclaxive values of inelastic rotations 
(actual '-'alue nultiplied by 6EI/L). 

Ro cation Capacity 

Tlic inelastic lengths are 

1^ = 0.227, I3 = 0.54, I4 = 0.273, I5 = 0.181 

Then relative values of rotation capacities (actual value 
i'QultipliGd by 6EI/I) are 

= I X 0.227 X 50 X I = H. 5 

= i X 0.273 X 50 X I = 10.2 

= f X 0.181 X 50 X I = 6.9 (6.29) 

Se cond Cycle of A na lysis 

The additional constr ants in second cycle of 
analysis are 

4M^ + 2M^ - 3M5 + 77.70A 4 8.5 

2M^ + 6M4 - 4M5 + 113.1 A 4 - 10.2 

3 M^ + 4M4 - 6M5 + 120. 9 A 4 ^*9 

4M^ + 2M4 - 3M5 + 77.0 A ^ 0.0 

2M^ + 6M^ - 4M5 + 113.1 A ^ 0.0 

3M^ + 2M4 - 3M5 + 77.70A!>0.0 


(6.30) 



65 


Maxinun value of Xand corj^esponding nonent distribution, 
subjGc b Lo conatr-int Eq.(6,22) and Eq.(6,30) is given 
in Table 6 .'\. 

T/iBLE B.i). Moment Distribution at Ultinato load 
for Example 6.2 , 


Seebion 

1 

2 

3 

4 

5 

M 

-17,67 

8.1 

16.3 

-22.0 

22.0 

A 

0 

0 

0 

7.66 

6.9 

®pi 

0 

0 

0 

10.7 

6.9 

load factor - 

2.32 





lloduction in linit load is 12,2 percent. 

Example 6,3 

It IS required to find out the ultir ate lo^d for the 
frame in Eig,6,5(a). The seetion properties for all the 
cribical sections are 

0<^ 1.18, B = 6.0, El = 2900, 11^=10.0 

Eirs fa Gycl o of Analysis 

V/ith no restriction on inelestic rotation at 
Q upiL ti, c al sec blons, the ultiiiatc load analysis problem is 
Maximise Z = A 
SubjocL to bho constraints 




ee 


M2 - 211^ + +20 X = 0 

Mg-.2IT^+MQ +20 A = 0 

*- + M5 + I.?g ^ 0 

+ M5 +MQ-''Ig-H^Q+20A = 0 

and -10 10 

i = 1,... 10 (6.31) 

'Plio oolut Lin or Lhis problen ic 

= -10.0 5 M 2 =• -6.67 ; = 10.0 ; 1.1^ = -10.0 

= -3.35 Mg = -6.67 5 My = 10.0 ; Mg = - 10.0 

Mq =• 10.0 $ = 10.0 5 A = 1 .833 (6.32) 

Jnel.alic "Rjoabijns 

The distribution of M , Ivl ond tl arc tobul' ted 

HJm LrJ X 

below in T''blo 6 . 5 . 

TABIjlil 6.5 Bcsidur>l Mooent Cclculatioiis for Ex 'uple- 6.3 


3c c tion 

1 

M 

Ul 


-10.0 

^ei 


- 4.2 



- 5.8 

Ru 01 dual nonont 
ru v/ith ad just - 
, ed sign 

- 5.8 


3 

4 

10.0 

- 10.0 

10.9 

-14.88 

-0.9 

4.88 

0.9 

4.80 
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The. poaiUvo dirootlon of the rosiduol nonenta are 
•jhov/ii 111 rijj, 6.5(d). The equilibnun equations in terns 
oj the rc&iducil nonent distribution are 


n-j + 2n^ + 2n^ + 2n^ 

02+203 f 

Og + 2nY 

^4 •" 1^5 


+2nQ+ng+n., 0 

- °8 


using tliGoe relations, residual rnonents 02 > 


= 0 
= 0 
= 0 
= 0 

(6.33) 
an(3 Og 


arc 


02 = - 6.68 

O 3 = 3*54 

rig = 1 .34 (6.34) 


■''isuning the last hinge to foro at Section 1 , the 
incLitabJC rot'etions at Sections 1,2,5 aiid 6 should be zero 
and lienoo , rip 02 j and Og arc dependent oonponents of 
residue'll inonciib system. Thus using Eq.(6,33) 



9rn.| 

_ 0 - 


^'"1 2 . 

2 ■ 

0i»3 ’ 

vl 

1 

11 

'CO 

- 2 . 

d’^ " ' 

3n3 - - > 

an, 

3ng 

3n, _ ^ 

31110 




.. .. < _ 4 . 

a^^4 


9”6 , 

S"8 ' ‘ 



0nit- 

and ~ = -1 

9^4 

9nic 

’ d^'j~ ^ ’ 

0"5 , 

S”8 

(6.35) 


Rotation calculation for this problon la given In Table 6 . 6 . 




TAELE 6.6 Rotation Calculation of Example 6,3 
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jjbcprcssions for residual norients in terns of actual nonent 
fit sections 1 j 3 , 4 , 7 , 8,9 and 10 whor^ hinges hare forned are 

+ 2.293 A 

rio = -- = 21^ ~ + 20.016 A 

~ + 5.95 A 

~ *'■ 8»116 A 

ni^ * - Pig = M 4 - 2Urj + Mg + 23.736 A 

ng = -2my - mg = 2My - Mg - 15.626 A 

Illy = -My + 5.313 A 
^8 = % 5.0 A 

nig = -Mg + 4 .785 A 

m^Q= -M^q + 4.17 A (6.36) 

Using bheoG expressions for residual nonents inelastic 
rota b LOIS arc 

0 ^ = .. 6 M^ - 2 OM 5 + 8 M 4 + 170.17 A 

= -5M^ - IOM3 H 8M4 - AMrj + 2Mq - + 159.1 A 

0^ = -4M^ - 4M3 + 6M4 - 16 My + aig - 2 M^q + 185. 6 A 

0Q = -4M^ - 4 M 5 + 4M^ - SMy + 6 Mg - Mg - M^q + 11 7. 9 A 

0g = -2M^ - 2M3 + + Mg - 2Mg ,+ 30.1 A 

0 ^Q= - 2 M^ ~ 2 M 5 + 2 M 4 - 2 My + Mg - 2 M^q + 47 .5 A 

(6.37) 
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Iho 0 v^^luco nrc all rclativo values of ix clastic 
3 ct Lion ( bht 'hboLute value = 0L/6i;i), 

riot bLon Cap city 

Incl >otic curvature for uoiicnt H bGtv.roen M and 

T/r . 



= { 


C>, 



(M - Mv) 



8.5) 


(6.38) 


1‘'jg.6.5(c) , the inelastic lerxgthsaro 
= 0.329 ; 1 , = 0.15 ; ly = 0,329 

■^S ^'‘^5 } Ig = 0.30 ; i^Q= 0.45 ( 6 . 39 ) 

Thun rx.iative values ot rot^ation capacities arc- (absolute 
value huutipli^d by 6EI/L) 



i y 0.329 X 26.8 X 1 ,5 X I 


= 10.2 


®p 4 “ 

OpO 13.95 
©p9 = 9.30 

©p10= ■'3.95 (6.40) 

Second Oyclo of Analysis 

The linit analysis problen vath rotation conpatibility 
consbrainb is 

Mincjiiiso ; Z = A 
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SubjGct to the constraints 

(a) Equilibrium equations and yield condition, Eq,(6.2l) 

(b) e^< 

and ^ ® 

where 0^ and rre given by Eq,(6.37) and Eq.(6,40) 
r-.spect3 vely. 

'J^ho solution of this problem is tabulated below 
111 'I iblo 6.7 , 

TiU3IiE 6.7 


Section 

1 

2 



5 

6 

7 

8 

9 

10 

M 

-3.61 



-10.0 

-6 .0 

-4.0 

10.0 

6.83 

8.39- 

6.8C 

6EI0^ 

L 

0.0 

0.0 

3.89 

4.65 

0.0 

0.0 

1 .38 

0.0 

0.0 

0.0 

i 

0.0 



4.65 

1 0.0 

0.0 

12.16 

jO.O 

0,0 

.1 

0.0 


Load factor 

= 1.54 








Reduction in load factor duo to limited ductility is 
16,3 percent. 
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The following table gives a comparision of the ultimate 
load calculated wibh and without limited ductility for a few 
typical structures. 
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CHAPTER VII 
CONCLUSIONS 
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1 . Proposed nethod using linear prograriDing can 
be viiry cffoctively used to arrive at a unique solution 
in plastic analysis of structures satisfying oquilibriun, 
yield and inechanism conditions, 

2. For analysis or design of fx’ancs larger than 
about 5 to 4 indetcrninecy , one has to always resort to 
conputor technique. The proposed pethod is ideally suited 
to coiaputi^r, 

3. Strain hardening effect is a very important 
factor to be tnhen into account for raxional evaluation 
of linit lords. The proposed nethod very conveniently 
Lakes this into account. 

]■ . Por materials with limited ductility, rotation 
c onpatibilitj’' has to be satisfied. This jlso is convenien- 
tly incorporated into proposed nethod, 

5. Prom the examples worked out, it is se^-n that fora 
structure with linitcd ductility, the ultimate load is 
reduced up to the uxtent of 15 - 20 ^ 0 , However, this exfont 
of reduction depends on available ductility. 
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